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Abstract. The existence of quasi-bi-Hamiltonian structures for the Kepler problem is 
studied. We first relate the superintegrability of the system with the existence of two 
complex functions endowed with very interesting Poisson bracket properties and then we 
prove the existence of a quasi-bi-Hamiltonian structure by making use of these two functions. 

The paper can be considered as divided in two parts. In the first part a quasi-bi-Hamiltonian 
structure is obtained by making use of polar coordinates and in the second part a new 
quasi-bi-Hamiltonian structure is obtained by making use of the separability of the system 
in parabolic coordinates. 
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1 Introduction 

1.1 A geometric introduction 

In differential geometric terms, the phase space M of a Hamiltonian system is the 2n-dimensional 
cotangent bundle M = T*Q of the n-dimensional configuration space Q. Cotangent bundles are 
manifolds endowed, in a natural or canonical way, with a symplectic structure wq; if {(%) | i = 
2,... ,n} are local coordinates in Q and {{qj,Pj)', j = 1,2,... ,n} the induced coordinates in 
T*Q, then ujq is given by 

loq = dqj A dpj, coq = —dOo, 6 q = Pjdqj, 

(we write all the indices as subscripts and summation convention on the repeated index is used). 
Given a differentiable function F in T*Q, F = F{q,p), the vector field Xp uniquely defined as 
the solution of the equation 

i{XF)ujo = dF 

is called the Hamiltonian vector field of the function F. In particular, given a Hamiltonian 
H = H{q,p), the dynamics is given by the corresponding Hamiltonian vector field Xh, that is, 
i{XH)oJo = dH. 

A vector field E G X{T*Q) is Hamiltonian if there is a function H such that E = Xp, i-e., 
i(E)u;o = dH, and locally-Hamiltonian when i(E)a;o is a closed 1-form. This is equivalent to 
E-invariance of wq, i.e., Cpojo = 0. 

A system of differential equations is called bi-Hamiltonian if it can be written in two different 
ways in Hamiltonian form. Suppose a manifold M equipped with two different symplectic 

*This paper is a contribution to the Special Issue on Analytical Mechanics and Differential Geometry in honour 
of Sergio Benenti. The full collection is available at http://www.emis.de/journals/SIGMA/Benenti.html 
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structures ujq and cui. A vector field F on T*Q is said to be a bi-Hamiltonian vector field if it 
is Hamiltonian with respect to both symplectic structures, i.e., 

i(F)a;o = dH^ and i(F)wi = dHi. 

The two functions, Hq and Hi, are integrals of motion for F. A weaker form of bi-Hamiltonian 
system is when the only symplectic form is the first one (cji is a closed but nonsymplectic 
2-form). 

We point out that an important example of bi-Hamiltonian system is the rational harmonic 
oscillator (non-central harmonic oscillator with rational ratio of frequencies) [8] 

H = ]^{pI+ pI) + {rri^x^ + n^y'^). 

A symplectic form determines a Poisson bivector A that satisfies the vanishing of the Schouten 
bracket [A, A] = 0 (this property is equivalent to the Jacobi identity); we note that there exist 
non-constant rank Poisson structure not related with symplectic structures. The compatibility 
condition between two different Poisson structures, Aq and Ai, means that the linear combina¬ 
tion A;, = Aq — AAi is a Poisson pencil, that is, it is a Poisson bivector for every value of A; 
therefore the corresponding bracket {•, -Ia = {•, - jo — A{-, -ji is a pencil of Poisson brackets. 

Bi-Hamiltonian systems are systems endowed with very interesting properties but, in general, 
it is quite difficult to find a bi-Hamiltonian formulation for a given Hamiltonian vector field, 
and for this reason it is useful to introduce the concept of quasi-bi-Hamiltonian system. 

Definition 1. A vector field W on a symplectic manifold {M,oj) is called quasi-Hamiltonian if 
there exists a (nowhere-vanishing) function y, such that yX is a Hamiltonian vector field 

yX e Xh{M). 

Thus i{yX)u! = dh for some function h. 

We call y an integrating factor of the quasi-Hamiltonian system, because it is an integrating 
factor for the 1-form i{X)uj, and we note that in this case the function h is a first integral of X. 
Note that this condition can alternatively be written as as i{X){yuj) = dh, but the point is that 
the 2-form yuo is not closed in the general case. 

The scarcity of bi-Hamiltonian systems leads to relax the condition for the vector field being 
Hamiltonian to a simpler situation of quasi-Hamiltonian with respect to the second symplectic 
structure. 

Definition 2. A Hamiltonian vector field W on a symplectic manifold (M, cj) is called quasi-bi- 
Hamiltonian if there exist another symplectic structure oji, and a nowhere-vanishing function y, 
such that yX is a Hamiltonian vector field with respect to wi. 

This concept was first introduced in [4] in the particular case of systems with two degrees of 
freedom and it was quickly extended in [23, 24] for a higher-dimensional systems. Some recent 
papers considering properties of this particular class of systems are [1, 2, 3, 4, 5, 6, 8, 14, 15, 23, 
24, 25, 36]. 

The nondegeneracy of the canonical form loq provides a vector bundle isomorphism cJq of 
T{T*Q) on T*(T*Q), inducing an identification of vector fields and 1-forms on the phase space. 
A consequence is that the pair (wq, wi) determines a (1,1) tensor field R defined as R = u)q~^ouji, 
that is. 


uJi{X,Y) =0Jii{RX,Y), yX,Y£X{T*Q). 
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Note that in the definition of R the 2-form oji is necessarily closed but it can be nonsymplec- 
tic. If r is bi-Hamiltonian with respect to (wo,wi) then R is P-invariant, that is, C^R = 0, 
where C denotes the Lie derivative (this means that the characteristic polynomial of R is an 
invariant for P, and consequently the coefficients of the polynomial are constants of motion). 
The Nijenhuis tensor N^i of the tensor field R is defined by 

Nr{X, Y) = R^[X, Y] + [RX, RY] - R[RX, Y] - R[X, RY]. 

It has been proved that if P is a Hamiltonian dynamical system of the type described above and 
such that (i) The tensor Nr of R vanishes, (ii) The tensor field R has n distinct eigenfunctions 
(that is, they are maximally distinct), then the eigenfunctions of R are in involution and the 
system is therefore completely integrable [12, 16, 17, 18, 21], It is important to note that the 
eigenvalues of R are constants of motion for P even in the case that the two properties (i) and (ii) 
are not satisfied (but then the eigenfunctions are not in involution). 

It is known that the Liouville formalism characterize the Hamiltonians that are integrable 
but it does not provide a method for obtaining the integrals of motion; therefore it has been 
necessary to elaborate different methods for obtaining constants of motion (Hamilton-Jacobi 
separability, Lax pairs formalism, Noether symmetries, Hidden symmetries, etc); the existence of 
a bi-Hamiltonian structure with the above two mentioned properties (Nijenhuis torsion condition 
and maximally distinct eigenvalues) can be considered as method to establish the Liouville 
integrability of a system; because of this, these two properties are frequently included in the 
definition of a bi-Hamiltonian structure. 

Most of systems admitting a bi-Hamiltonian structure are separable systems; so these two 
properties (separability and double Hamiltonian structure) are properties very close related 
(see [18] for a detailed discussion of this question and [35] for the case of multiple separability). 
Quasi-bi-Hamiltonian systems are very less known than the bi-Hamiltonian ones but it seems 
that they are also related with separability. Nevertheless, in this case the tensor field R is not 
P-invariant and the eigenvalues of the tensor field R are not constants of motion. An interesting 
property is that it was shown in [4] that in the particular case of two degrees of freedom the 
function / det R is a constant of the motion. 

The potential of the Kepler problem is spherically symmetric and therefore it admits al¬ 
ternative Lagrangians (the existence of alternative Lagrangians for central potentials is stud¬ 
ied in [13, 20, 26]); recall that if there exist alternative Lagrangian descriptions then one can 
find non-Noether constants of motion [7]. This system has been studied as a bi-Hamilltonian 
system by making use of different approaches; Rauch-Wojciechowski proved the existence of 
a bi-Hamiltonian formulation but introducing an extra variable so that the phase space is odd¬ 
dimensional and the Poisson brackets are degenerate [34] and more recently [19] a bi-Hamiltonian 
formulation for the perturbed Kepler problem has also been studied by making use of Delaunay- 
type variables. 

1.2 Structure and purpose of the paper 

Now in this paper we will analyze certain properties of the Kepler problem related with the 
existence of quasi-bi-Hamiltonian structures. The following two points summarize the contents 
of the paper. 

• First, we will study the existence of certain complex functions with interesting Poisson 
properties and then we will prove that the superintegrability of the system is directly 
related with the properties of these complex functions. As is well-known this system is 
multi-separable (it separates in both polar and parabolic coordinates); so we first present 
the study making use of polar coordinates and then we undertake a similar study by 
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employing parabolic coordinates (the parabolic complex functions are different from the 
polar ones). 

• Second, we prove that the above mentioned complex functions determine the existence 
of several quasi-bi-Hamiltonian structures. This is done in two steps: first with complex 
2-forms and then with several real 2-forms. The properties of these geometric structures 
and of the associated recursion operators are analyzed. 

It is important to note that this study is concerned with the existence of quasi-bi-Hamil¬ 
tonian structures (instead of bi-Hamiltonian). So we recall that if z(F)u;i = XdHi then 
Xyuji = dX a dHi 7 ^ 0. Consequently the tensor field R is not F-invariant and the 
eigenvalues of R are not constants of motion. 

We must clearly say that the structures obtained by this method (wedge product of the 
differentials of complex functions) do not satisfy the above mentioned Nijenhuis torsion con¬ 
dition. So perhaps it is convenient to name them as weak quasi-bi-Hamiltonian structures (in 
opposition to strong structures satisfying the Nijenhuis condition). Nevertheless, the purpose 
in this paper is not to prove the integrability of a system as consequence of a bi-Hamiltonian 
structure, since it is perfectly known that the Kepler problem is not only integrable but also 
superintegrable. The purpose is to study new and interesting properties of the Kepler problem. 
In fact, it has been proved that if a dynamical vector field satisfies certain properties (existence 
of canonoid transformations [9, 10] or existence of non-symplectic symmetries [27, 28]) then it is 
Hamiltonian with respect to two different structures without satisfying necessarily the Nijenhuis 
condition. 

The structure of the paper is as follows: In Section 2 we study the Kepler problem by making 
use of polar coordinates (r, </)); we relate the superintegrability of the system with the existence 
of two complex functions Mr and endowed with very interesting Poisson bracket properties 
and then we prove the existence of a quasi-bi-Hamiltonian structure making use of these two 
functions. Then in Section 3 we consider once more the same system but in terms of parabolic 
coordinates (a, b) and we obtain a new quasi-bi-Hamiltonian structure (different to the previous 
one) making use of a similar technique but with new complex functions Ma and Mi,. Finally in 
Section 4 we make some final comments. 


2 Prom superintegrability to quasi-bi-Hamiltonian structures 

After these rather general comments we restrict our study to the Kepler problem in the Euclidean 
plane that, as it is well known, is superintegrable and multiseparable (polar and parabolic 
coordinates). 

Let us first notice that in some cases the two-dimensional Euclidean systems possess certain 
interesting properties. Eor example, if the potential V{x,y) takes the form V = A{u) + B{v), 
u = x + y,v = x — y, then it admits a new Hamiltonian structure (and also a Lax pair) [22]; 
unfortunately the new structure is in most of cases constant (we mean that the new Poisson 
bracket is determined by a symplectic form as loi = dx A dpy + dy A dpx). 

In the next paragraphs we will relate the existence of bi-Hamiltonian structures with the 
properties of the two complex functions with interesting Poisson bracket properties. 


2.1 Complex functions and snperintegrabilty 

It is well known that the Hamiltonian of the two-dimensional Kepler problem 


Hk = 


1 

2 



Vk = 


9 


+ Vk 


r 


0 < g eR, 
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is Hamilton-Jacobi (H-J) separable in polar coordinates (r, (/>) and it is, therefore, Liouville 
integrable with the angular momentum J 2 = P(p as the associate constant of motion. Moreover, 
it is also known that it is a super-integrable system with the two components of the Laplace- 
Runge-Lenz vector as additional integrals of motion. Now we will prove that this property of 
superintegrability can be related to the existence of certain complex functions with interesting 
Poisson bracket properties. 

Let us denote by M^j and j = 1,2, the following real functions 

p2 

Mri = PrP 4 >, Mr 2 = g - and = cos(/>, N^2 = sin cj). 

Then we have the following properties 

(i) ^Mrl = {Mrl, Hk} = -XMr2, ^Mr2 = {Mr2, Hk} = AM,i, 

at at 

(ii) = {Xlcpi, Hk} =—XNfj,2, = {Hfj,2, Hk} = XN^i, 

where A denote the following function 

The property (ii), representing the behaviour of the angular functions is true for all the 
central potentials H(r); but the property (i), behaviour of the functions Mrj, is specific of the 
potential of the Kepler problem. 

Consider next the complexification of the linear space of functions on the manifold and extend 
by bilinearity the Poisson bracket. If we denote and the complex functions 

Mr = Mrl + iMr2 , = Nfpi + iN^2 , 

then they have the following properties 

{Mr, Hk} = iXMr, Hk} = iXNip, 

and consequently the Poisson bracket of MrN^ with the Kepler Hamiltonian vanishes 
[MrN;, Hk} = {Mr, Hk}N; + Hk} = (iAM,)iV; + Mr{-iXN;) = 0. 

We can summarize this result in the following proposition. 

Proposition 1. Let us consider the Hamiltonian of the Kepler problem 

Hk = ^(^Pr + ^^ +Vk, Vk = -^, 

Then, the complex function J34 defined as 
J34 = MrN; 

is a (complex) constant of the motion. 

Of course J34 determines two real first-integrals 


J34 = J3 + iJ4:, {J3 ,Hk}=0, {J4,Hk}=0, 
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whose coordinate expressions turn out to be 

J 3 = Re( J 34 ) = PrPrj) COS (j) - - sin (/) + g sin cj), 

p2 

J 4 = Im( J 34 ) = PrPtj, sin (j) -\—- cos (f) — g cos cj). 


That is, the two functions J 3 and J 4 are just the two components of the two-dimensional Laplace- 
Runge-Lenz vector. 

Summarizing, we have got two interesting properties. First, the superintegrability of the 
Kepler problem is directly related with the existence of two complex functions whose Poisson 
brackets with the Hamiltonian are proportional with a common complex factor to themselves, 
and second, the two components of the Laplace-Runge-Lenz vector appear as the real and 
imaginary parts of the complex first-integral of motion. Remark that is a complex function 
of constant modulus one, while the modulus of is a polynomial of degree four in the momenta 
given by 


MrU; 


{PrP<t>f + 



2plHK + g^. 


2.2 Complex functions and quasi-bi-Hamiltonian structures 

Let us denote by I 34 the (complex) Hamiltonian vector field of J 34 

^(^34)1^0 = dJ 34 , 


that obviously satisfies Y^ 4 {Hk) = {Hr-, Jm} = 0 , and by and Y^ the Hamiltonian vector 
fields of Mr and Nfj,: 

‘i'(Yfp')(jjQ — dN(j). 


Their local coordinate expressions are, respectively, given by 


d 


.P4> 


-P<t>Q^ + \Pr ^ ) PiA, ^2 ^ 


d 


d 


d(j) 


dpr 


and 

d 

Ya = (sin 6 — 1 cos 6) -—. 

dp4> 

Then, the vector field T 34 appears as a linear combination of and Y ^; more specifically we 
have 


Y34 = N;Yr + MrY;=Y + Y', Y = N;Yr, Y' = MrY;. 

The vector field ^34 is certainly a symmetry of the Hamiltonian system {T*Q,ujq, Hk), but the 
two vector fields, Y and Y\ are neither symmetries of the symplectic form uoq (that is, XyojQ 7 ^ 0 
and Xy'UJq 7 ^ 0) nor symmetries of the Hamiltonian (that is, XyPix 7 ^ 0 and Xy/Hx 7 ^ 0). 
Moreover, remark that they are not symmetries of the dynamics, because 

[T, F,^] 7^ 0, [y', F;^] 7^ 0, i{rK)u:o = dHx- 


More specifically: 
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Proposition 2. The Lie bracket of the dynamical vector field Tk with Y is given by 

[rK,Y] = iJsax, 

where X\ is the Hamiltonian vector field of the function A. 

Proof. A direct computation leads to 


[r^, K] = TK{N;)Yr + n;[Tk, k.] = -iXN^Yr + n; • 

where we have used that the Lie bracket of two Hamiltonian vector fields satisfies [Xf,Xg] = 
—Xygy Note also that the Hamiltonian vector field of a product fg is given hy Xfg = 
fXg + gXf, and then the above Lie bracket becomes 

[Fk, Y] = -iXN^Yr + iN;{XxM^) = -iXN^Yr + iN;{XYr + M,Aa) = i{MrN;)Xx. ■ 

The vector field Xx on the right hand side represents an obstruction for Y to be a dynamical 
symmetry. Only when A be a numerical constant the vector field Y (and also Y') is a dynamical 
symmetry of 

In the following H will denote the complex 2-form defined as 
Q = dMrAdN^. 

The two complex 2-forms wy and Uy obtained by Lie derivative of ujq, i.e., 

TyCUO = Ty'COq = UJyi 

are such 

CyoJo = iyidujo) + d(iycJo) = d{iyuJo) = d{N^dMr) = —Ft, 

CyiUjQ = iyiifilOo) -\- dfiyiUlfi) = d{iyiU)fi) = d^Mj-dN^) = H. 

Using the preceding results we can prove: 

Proposition 3. The Hamiltonian vector field Tx of the Kepler problem is a quasi-Hamiltonian 
system with respect to the complex 2-form H. 

Proof. The contraction of the vector field Tk with the complex 2-form H gives 
i{rK)n = rK{Mr)dN; - ri^(iv;)dM„ 
and recalling that 

TKiMr) = {Mr,HK} = iAM„ TKiN;) = {N;,Hk} = -iXX;, 

we arrive to 

i{rK)Fl = {iXMr)dNl + {iXN^)dMr = iXd{MrN^). ■ 

The complex 2-form H can be written as 
H = -\- iFl2 

where the two real 2-forms, Hi = Re(H) and H 2 = Im(H), take the form 

Hi = dMri A dNipi dMr 2 A dNrp 2 = d{prP^) A d{cos (/)) -\- d (g — —^ A (i(sin (f) 
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^}2 = —dMri A dN^2 + dMr2 A dN^i = —d{prP^) A d{sm(f)) + d ^51- - j A d{cos(j)) 

= Pudr Ad(j) + I32zd4> A dpr + hidcf) A dp^ 
with aij and fdij being given by 

pI 


ai2 = ( ^ I cos (/), 023 = sin (/), 024 = Pr sin 0 + 2 cos (/>, 


and 


2 

/3i 2 = - f ^ j sin (/), (323 = Pcpcoscj), /?24 = Pr cos (/> - 2 sin (/>. 


Then we have 

i(F/f)r2i = —XdJi, i(F^)r22 = XdJs, 


what means that F^^ is also quasi-bi-Hamiltonian with respect to the two real 2-forms (a;o,Hi) 
or ( 020 ,^ 2 ). 

Remark that the complex 2-form H is well defined but it is not symplectic. In fact, from 
the above expressions in coordinates we have Hi A Hi = 0, H 2 A H 2 = 0, and Hi A H 2 = 0, and 
therefore we obtain 


H A H — (Hi A Hi — H 2 A H 2 ) -I- 2iHi A H 2 — 0. 


The distribution defined by the kernel of H, that is two-dimensional, is given by 
KerH = {/iZi + / 2 Z 2 |/i,/ 2 : xM^^C}, 

where the vector fields Zi and Z 2 are 

d . d . d 

Z\ = (023 + i(d23)-^ + {oii2 + i(3i2)-^-, Z2 = (024 + + (012 + i(di2) 

Therefore it satisfies 


d 

dp<i>' 


[KerH,Fi^] C KerH. 

That is, Fi^ preserves the distribution KerH. 

If 13 and F 4 are the Hamiltonian vector fields (with respect to the canonical symplectic 
form cuo) of the first integrals J 3 and J 4 , then the dynamical vector field Fx is orthogonal to F 4 
with respect to the structure Hi and it is also orthogonal to Ts with respect to the structure H 2 , 
that is. 


i{TK)i{Yi)n^ = 0, i{TK)i{Y^)^2 = 0. 


Just to close the section we remark that had we applied this technique to the isotropic two- 
dimensional harmonic oscillator with frequency a we had obtained the function as 


Mr = ( ) + i 


{ N 

2 ^ 22 
Pr — ~2 + CU r 


(the angular function would be the same) and the constants so obtained are but the compo¬ 
nents of the Fradkin tensor. This shows that the harmonic oscillator is an example of dynamical 
system both bi-Hamiltonian and quasi-bi-Hamiltonian. 
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2.3 Recursion operators and some comments 

The bi-Hamiltonian structure determines a complex recursion operator R defined as 

Q{X,Y) =ujo{RX,Y), yX,YeX{T*Q). 

But as H and R are complex, we can introduce two real recursion operator Ri and R 2 defined 
as 


Hi(x, y) = cooiRiX, y), Q2{x, y) = ioo{R2X, y). 

We recall that ojq is the map wq- 3i{T*Q) —)■ A^{T*Q) given by contraction, that is = 

i{X)uJo, and then the nondegenerate character of ujq means that the map wq is a bijection. Using 
this notation we can write the two operators Ri and R 2 as follows 

Rl = UJQ ^ o Hi, R2 = OJQ ^ o H2. 


Then we have the following properties 

(i) The coordinates expressions of Ri and R 2 are 


d 


d d d 

023^ + “24^ + ai2x— 
C/(p 


Rl = -0127,— (8) dr -h I .. 

op^ I or 

d d 

+ 023 X- ® dpr -h 024 7,- <8 dps 

op^ op^ 


® d<p 


and 


d 

R 2 = -/ 5 i 2 ^— ®dr + 

op<t> 


d d d 

hzTT- + - 

or ocp OPr 


® d(j) 


ad O d 

P 23 ® dpr -|- P24 ^ ® dpfj). 

OP<f> OPcP 


(ii) Rl and R 2 have two different eigenvalues doubly degenerate and one of them is null (that 
is, Ai = A 2 = 0, A 3 = A 4 / 0). Therefore we have 

det[i?i] = det[ii 2 ] = 0 , 

what is a consequence of the singular character of Hi and H 2 . 

We close this section summarizing the situation we have arrived. We have first introduced 
two complex functions, and iV,^, mainly because of the behaviour of their Poisson brackets. 
Then we have proved that they are interesting for two reasons; first because they determine 
the existence of superintegrability (existence of additional constants of motion) and second 
because they determine quasi-bi-Hamiltonian structures (first complex (wo,H) and then real 
(cuo, Hi, H2)). 

Concerning the first point, in this case the additional constants of motion are just the com¬ 
ponents of the Runge-Lenz vector (that have been highly studied making use of different ap¬ 
proaches). Now we have arrived to a new property; they can also be obtained as a consequence 
of this complex formalism. 

Concerning the second point, the two complex functions and determine the above 
mentioned geometric structures (first complex and then real) but unfortunately they are dege¬ 
nerated (we recall that Hi A Hi = 0, H 2 A H 2 = 0). This can be considered as a limitation of 
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these geometric structures. If a bi-Hamiltonians structure satisfies all the appropriate proper¬ 
ties (that is, symplectic forms, vanishing of the Nijenhuis torsion of the recursion operator R, 
diagonalizable recursion operator R with functionally independent real eigenvalues) then it de¬ 
termines the Liouville integrability of the system. In fact, the aim of the approach presented 
in this paper is not to prove the integrability of a system as consequence of a bi-Hamiltonian 
structure as we start with a system known to be not only integrable but also superintegrable. 
The existence of (111,112) must be considered, not as a method for arriving to the integrability 
of the system, but as a new and interesting property of the Kepler problem (for the moment 
only of the two-dimensional system, the generalization to the three-dimensional case must be 
considered as an open question). 


3 New complex functions and new quasi-bi-Hamiltonian 
structures 

In this section we will study the existence of new bi-Hamiltonian structures for the Kepler 
dynamics by making use of parabolic coordinates (a, b) defined as 

X = 2 , y = oh. 

Of course all previous results can be translated to this new language in such a way that the 
functions M^i and Mr 2 are now given by 

{apb - bpa){apa + hpi,) {oph-bpaf 

Mrl = - / -, Mr2 = - , - g, 


\/a2~+~P ’ 

while functions N^i and N ^2 become 




a? - 6 ^ 


N^2 = 


2ab 


+ b‘^ 


But the important point is that the behaviour of the Kepler Hamiltonian in these coordinates 
will permit us to obtain new results different from the previous ones. 


3.1 Complex functions and superintegrability 

The general form of a natural Euclidean Hamiltonian is 

1 


H = 


pI+pI 


+ V{a,b). 


2m \a‘^ + b‘^ 

in such a way that if the potential V is of the form 
A{a) + B{b) 


Via,b) = 


a? + 


then the Hamiltonian is Hamilton-Jacobi separable and it is, therefore, Liouville integrable with 
the following quadratic function 

1 f a? B — b^A\ 

•^2 = ( q 2 ^ 52) - bPa){opb + bpa) + 2 ^ J 


as the second constant of motion (the first one is the Hamiltonian itself). 
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For simplifying the following expressions we introduce the following notation: 


J = apb- bpa, Px = 


gpa - bpb 
a'^ + b'^ ’ 


apb + bpa 
a? + 6^ 


The Hamiltonian of the Kepler problem when written in parabolic coordinates is 


Hk = 


1 (pI+pI 


2 I a 2 + 62 


+ Vk, Vk = - 


:2 + 62 ’ 


and the Kepler dynamics is given by the following vector field 


rx = 


Pa \ d 


d , f pl+pl-‘ig\ d fpl + pl-2g\ d 


a^ + b‘^J da + b"^ J db \ (a2-|-62)2 j Qp^ 


- 1 " 


(a2 -|-62)2 J Qp^'’ 


in such a way that, as loq = da A dpa + db A dpb, we have 
i{TK){da A dpa + db A dpb) = dHx- 

The Hamiltonian Hk is Hamilton-Jacobi separable in coordinates (a, 6) and the associated 
quadratic constant of motion is the component Rx of the Laplace-Runge-Lenz vector 


Rx — J Py 9 


a? ) ' 


Let us now introduce the functions Maj and Mbj, j = 1,2, defined by 


Mai = 




Ma2 — 


2 ga - Jpb 

v/FTF’ 


Mbi = 


r 2 + 62 ' 


Mb2 = 


2gb + J Pa 

a/F+F 


Then, the important property is that the Poisson bracket of the function Mai with Hk is 
proportional to Ma 2 while the Poisson bracket of Ma 2 with Hk is proportional to Mai, but with 
the opposite sign; 

{Mal,HK] = —\Ma2, {Ma2, Hk] = XMal, 

and the same is true for the functions Mbi and Mb 2 , 

{Mbi,HK] = —XMb2, {Mb2,HK} = XMbi, 

where now A denotes the following function 

apb- bpa 

(a 2 + 62)2 • 

Therefore the two complex functions Ma and Mb defined as 


Ma = Mai + iMa2, Mb = Mbi + iMb2, 


satisfy 


{Ma, Hk} = iXMa, {Mb, Hk} = iXMb- 
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Proposition 4 . The complex function defined as 
K34 = MaM; 

is a {complex) constant of the motion for the dynamics of the Kepler problem described by the 
Hamiltonian (1). 

We omit the proof because it is quite similar to the proof of the previous Proposition 1. 
Note that the modulus of the complex functions Ma and are given by 

MaM: = 2 {J^Hk - gR. + g^), MbMl = 2 {J^Hk + gR. + g ^), 

and then 


MaM: + M^M: = + g^) ■ 

Of course the complex function determines two real functions that are first integrals for 
the Kepler problem 

^34 = K 3 + ii^4, {i^3, Hk) = 0 , {Ki, Hk} = 0 , 

with K 3 and iF 4 given by 

Ks = Re(K34) = MaiMbi + = JPx+g > 

-K 4 = Im(R'34) = Ma2Mbi — MalMiy2 = —2J^Hk- 

Remark that the function K 3 is the other Laplace-Runge-Lenz constant, while -K 4 , that is 
a fourth order polynomial in the momenta, determines the angular momentum J as a factor. 


3.2 Complex functions and quasi-bi-Hamiltonian structnres 

First we recall that the complex functions Ma y are given by 


Ma 


Jpa \ ^ ■ ( ‘2^ga - Jpb '] 

y/affiRfi J V \/a2 + 52 J 


Mb 


(M£LJ]+i(^4±2ts.) 


Let us now denote by Z 34 the Hamiltonian vector field of the function K 34 , i.e., i{Z 34 )u!Q = 
dZ 34 , such that Z 34 {Hk) = 0, and by Za and Zb the Hamiltonian vector fields of the complex 
functions Ma and Mb, that is, 


i{Za)ujQ = dMa, i{Zb)u}o = dMb- 


Their coordinate expressions are given by 


= 


dMa\ d 

' 7W + 


dMa\ d 


dMa \ 
da J dpa 


^ - 7 

r\i I CA ^<l\ \ ^^Obl'i 

ob J opb 


dpa J da \ dpb J db 
with Zai and Za 2 given by 

V 1 fr ni. \ ^ , 9 {apa + hpb) f^d d 

^al = / ^ ( {apb - + ^'P°‘Rh ^2 . m 1 - °^Pa — 


Za2 — 


1 

\/a 2 ~+~P 


Mb 


+ 62 


dpa 


L r\ 

dpb 


1.9 d {apa + bpb)pb -2gbf d d 

6pfe^ + {bpa - 2apb)-^ - o . ,n - I -b- -h a- 


db 


+ 62 


dpa dpb 
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and 


^6 = 


dMb\ d 

7^ + 


dpa J da 
with Zbi and Zb 2 given by 
1 


dMb\ d fdMb\ d 


dpb J db \ da J dp, 


dMb\ d _ 

— — ^bl ~r l^b2i 


db J dpb 


Zbi — 

Zb2 = 


1 


, <9 , .d {apa + bpb) 

-bp,- + (2ap,-bp„)-- 


, d d 

bpb^ - 

^Pa ^Pb 


! nb ^ ^ + bPb)Pa - 2£/a 

(“» - -- 


b-- — 

dpa dpb 


\/aF+~P 

Now recalling that 

dZu = d{MaMn = M^d{Ma) + Mad{M*b), 

we obtain 

Z 34 = M*bZa + MaZ*b =Z + Z', where 
The following proposition is similar to that of Proposition 2 and we omit the proof. 


^ = M^Za, 


Z' = MaZ*b. 


Proposition 5. The Lie bracket of the Kepler dynamical vector field Tk with the vector field Z 
is given by 

Z\ = iK^iXx, 

where X\ is the Hamiltonian vector field of X solution of the equation i{Xx)uJo = dX. 

In the following we will denote by Llab the complex 2-form defined as 


flab = dMa A dM^ 
JPa 


= d 


_ V\FF~+~P 


i 


2ga - Jpb 


A d 


JPb 




— i 


2gb + J Pa 

\4?T~P 


Then the two 2-forms ojz and obtained by Lie derivation of wq with respect to Z and Z' are 
given by 

T-zWlO = LOz = — flab, Tz'WJQ = Lo'z = flab- 

Proposition 6. The Hamiltonian vector field Tk of the Kepler problem is quasi-Hamiltonian 
with respect to the complex 2-form flab- 

Proof. This can be proved by a direct computation 

iiTK)flab = TK{Ma)dM^ - TK{M^)dMa = {iXMa)dM*b + {iXM^)dMa = iXd{MaM^). m 

The complex 2-form flab can be decomposed as 

flab — flabl T iflab2, 


where the two real 2-forms, flabi = Re(Dafe) and flab 2 = Iiii(Dafe); take the form 

^abl = dMal A dMbl -|- dMa2 A dMb2, flab2 = —dMal A dMb2 + dMa2 A dMbl, 
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Considering the real and imaginary parts we obtain 
i(F/f)na61 = —XdKi, iiJ'K)^ah2 = ^dK^, 

what means that F^ is also quasi-bi-Hamiltonian with respect to the two real 2-forms {ujq, ^abi) 
and (a;o,H„fe 2 )- 

Once more we obtain that the factor A determines that the system is quasi-bi-Hamiltonian 
instead of just bi-Hamiltonian. 

The complex 2-form ilab is well defined but it is not symplectic. The kernel is two-dimensional 
and it is invariant under the action of Fx 


[KerHafejF/^] C Kerflaj,. 


The coordinate expressions of Oabi and Qab 2 are 
2 J 


^abl — 


^ab 2 — 


(a 2 + 62)2 

(a2+62)2 


{aisda A dpa + ai^da A dpb -|- a23db A dpa -|- a2Adb A dpb + a^^dpa A dpb), 
(Pisda A dpa + fiiida A dpb + l323db A dpa + ^Adb A dpb), 


with aij and fdij being given by 

ai 3 = [2gb - apaPb - bpl)b, 
023 = + apaPb + bpl)a, 

034 = 2 J -|- 6^), 


ai 4 = -{ 2 ga - apl - bpaPb)b, 
024 = ( 2 ga - apl - bpaPb)a, 


^13 = {2abpa - a^pb - b‘^Pb)b, (du = {2abpb - o^pa - b‘^Pa)b, 

1^23 = {-2abpa + a^pb + b‘^Pb)a, (^24 = {-2abpb - o^pa - b^Pa)a. 


We close this section with the following properties: 

(i) The two real 2-forms are not symplectic. In fact we have verified that Hi A fli = 0, 
H 2 A H 2 = 0, and also Hi A H 2 = 0. 

(ii) These two 2-forms, Habi and ^ab 2 i determine two recursion operators ((1,1) tensor fields) 
Rabi and Rab 2 defined as 


H„fei(W,y) =a;o(RawX,y), = u^{R,b 2 X,Y), 


or in an equivalent way 


Rabl — ® H(jfti, Rab2 — ^0 ® ^ab2- 


-1 


As in Section 2.3, a consequence of the singular character of Habi and ^ab 2 is that 
det[Rabi] = det[Ra62] = 0. 


(hi) If we denote by Z 3 and ^4 the Hamiltonian vector fields (with respect to the canonical 
symplectic form wq) of the integrals K 3 and A' 4 , then the dynamical vector field F/^ is 
orthogonal to Z 4 with respect to the structure Hi and it is also orthogonal to Z 3 with 
respect to the structure H 2 , that is. 


*(rx)*(h4)Ha6i = 0, ax K)'iiY3)Q.ab2 = 0. 
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4 Final comments 

The Kepler problem is separable in two different coordinate systems, and because of this, it 
is superintegrable with quadratic integrals of motion. Now we have proved that this super- 
integrability is directly related with the existence of certain complex functions possessing very 
interesting Poisson bracket properties and also that these functions are also related with the 
existence of quasi-bi-Hamiltonian structures. 

We finalize with some questions for future work. First, as stated in the Introduction, quasi- 
bi-Hamiltonian structures is a matter that still remain as slightly studied (in contrast to the 
bi-Hamiltonian systems); so the particular case of the Kepler problem can be a good motivation 
to undertake a better study of these systems. Second, the existence of these complex functions 
is not a specific characteristic of the Kepler problem; in fact, it has been proved that the 
superintegrability of certain systems recently studied (as the isotonic oscillator or the TTW 
or PW systems) [11, 29, 30, 31, 32, 33] is also related with such a class of complex functions. 
Therefore, an interesting open question is whether these other superintegrable systems are also 
endowed with bi-Hamiltonian structures or with quasi-bi-Hamiltonian structures; the starting 
point for this study must be a deeper analysis of the properties of these complex functions 
making use of the geometric (symplectic) formalism as an approach. 
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